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The shape transformations of fluid membranes induced by curved protein rods are studied using 
meshless membrane simulations. The rod assembly at low rod density induces a flat membrane 
tube and oblate vesicle. It is found that the polyhedral shapes are stabilized at high rod densities. 

The discrete shape transition between triangular and buckled discoidal tubes is obtained and their 
curvature energies are analyzed by a simple geometric model. For vesicles, triangular hosohedron 
and elliptic-disk shapes are formed in equilibrium, whereas tetrahedral and triangular prism shapes 
are obtained as metastable states. 


I. INTRODUCTION 

Single-component fluid vesicles exhibit various mor¬ 
phologies such as stomatocyte, pear, and starfish 
shapes. These shapes can be understood in terms of 
the membrane bending free energy with area-difference 
elasticityji^— A red-blood-cell shape, discocyte, can be 
obtained by the energy minimization of the bending en¬ 
ergy with fixed area and volume. However, many cell 
organella have much more complicated shapes depend¬ 
ing on their functions, and their shapes cannot be re¬ 
produced by single-component lipid membranes. For 
example, the Golgi apparatus consists of stacks of dis¬ 
coidal membranes. These membrane shapes are regu¬ 
lated by various proteinsHowever, the mechanism 
of the shape regulation has not been clarified. 

In particular, protein domains called BAR 
(Bin/Amphiphysin/Rvs) domainsi^^— have recently 
been receiving increased attention experimentally and 
theoretically. The BAR domains consist of a banana¬ 
shaped dimer and are categorized into subsets of unique 
families including “classical” BARs, F-BARs (Fes/CIP4 
homology-BAR), and I-BARs (Inverse-BAR). The BAR 
domain absorbs onto the biomembrane and changes the 
local curvature. The extension of membrane tubes from 
liposomes and the formation of a cylindrical scaffold have 
been experimentally observedJ^^— Dysfunctions 
of the BAR proteins are considered to be implicated 
in neurodegenerative, cardiovascular, and neoplastic 
diseases. 

Membrane tabulation and budding can be generated 
by the adhesion of objects inducing isotropic sponta¬ 
neous curvature such as spherical colloids and coni¬ 
cal transmembrane moleculesi^ In the last decade, the 
membrane-mediated interactions between the isotropic 
adhesives were intensively investigated<^“— In contrast, 
the BAR domains are banana-shaped and generate an 
anisotropic curvature. This anisotropic spontaneous cur¬ 
vature is considered to give pronounced effects to con¬ 
trol local membrane shapes. Recently, the classical 
Canham-Helfrich curvature free energy^ii was extended 
to anisotropic curvaturesThe adsorption and as¬ 
sembly of the BAR domains were investigated using 


atomic and coarse-grained molecular simulations 
Formation of discoidal vesicles and tubulation were 
simulated using a dynamically triangulated membrane 
mode l^^i^^ and meshless membrane modelsi^^— Despite 
these numerous advancements, the physics of membrane 
shape deformation caused by anisotropic curvature is not 
completely understood. 

In this paper, we describe shape transformations of 
membrane tubes and vesicles using an implicit-solvent 
meshless membrane model.— In particular, we focus 
on the membrane shapes at high protein densities. A 
BAR domain is modeled as a banana-shaped rod, and it 
is assumed to be strongly adsorbed onto the membrane. 
In order to investigate the membrane-curvature-mediated 
interactions, no direct attractive interaction is considered 
between the rods. The membrane-mediated interactions 
induce the assembly of the protein rodsj^ We will show 
the assembly of protein rods induces the formation of 
polygonal tubes and polyhedral vesicles. Polyhedral vesi¬ 
cles were previously reported for gel-phase membranes^, 
phase-separated membranes^^— , and fluid membranes 
with the accumulation of specific lipids or defects on the 
polyhedral edges or vertices— Here, we demonstrate 
that the polyhedral vesicles can be stabilized by the pro¬ 
tein rods. 

In Sec. m the simulation model and method are pro¬ 
vided. In Sec. nni the shape transformations of the mem¬ 
brane tubes induced by the protein rods are described. In 
Sec. IIVI a simple geometric model is proposed to under¬ 
stand the results of Sec. IHII In Sec. El the shape trans¬ 
formations of vesicles are described. The summary and 
discussion are given in Sec. I VII 


II. SIMULATION MODEL AND METHOD 

The details of the meshless membrane model and pro¬ 
tein rods are described in Ref. 113 and Ref. [sl, respec¬ 
tively. Here, we briefly describe the model. A fluid mem¬ 
brane is represented by a self-assembled one-layer sheet 
of N particles. The position and orientational vectors of 
the f-th particle are and Ui, respectively. The mem¬ 
brane particles interact with each other via a potential 
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U — t^rep“t“C^att • TllG potential U^ep IS an GX 

cludGd volumG interaction with a diameter cr for all pairs 
of particles. Solvent is implicitly accounted for by an 
effective attractive potential 17att • It is a pairwise attrac¬ 
tive potential at low particle density while the attraction 
is smoothly truncated at high density. This truncation 
allows formation of a fluid membrane in a wide range of 
parameters. The bending and tilt potentials are given by 

Ifbend — (^bend/2) and 

Utat = (fctiit/2)E*<j[(u* ■ + (uj ■ rij)^]wcAn,j), 

respectively, where nj = - r^, Vij = \Yi^j\, Vij = 

Yijjrij, and Wevi^ij) is a weight function. The sponta¬ 
neous curvature Co of the membrane is given by Cocr = 
Cbd/24S In this study, Cq = 0 except for the membrane 
particles belonging to the protein rods. 

The protein rod is modeled as a linear chain of Ngg 
membrane particles. We use Ngg = 10, which corre¬ 
sponds to the typical aspect ratio of the BAR domains. 
The BAR domain width is approximately 2 nm, and 
the length ranges from 13 to 27 nmj^^ The protein rod 
has spontaneous curvature Crod along the rod axis and 
zero spontaneous curvature perpendicular to the rod axis. 
Hereafter, we call the membrane particle consisting of 
the protein rod a protein particle. The protein particles 
in each protein rod are connected by a bond potential 
ffrbond/^B^ — (^rbond/2o' )(Ti-\-l^i /rod) ■ The bend¬ 
ing potential is given by t/rbend/feeT’ = (fcrbend/2)(fi+i,i • 
— Cr)^, where Cr = 1 — (Crod/rod)^/2. For bonded 
pairs of protein particles, the corresponding spontaneous 
curvature and four times larger values of the coefficients 
of t/bend and Utat are employed to prevent the rod from 
bending tangentially to the membrane. 

We employ the parameter sets used in Ref. The 
membrane has mechanical properties that are typical 
for lipid membranes: the bending rigidity n/kUT = 
15 ± 1, the area of the tensionless membrane per particle 
= 1.2778 ± 0.0002, the area compression modu¬ 
lus Ka(t^/ kUT = 83.1 ± 0.4, and the edge line tension 
Ta/ksT = 5.73 ± 0.04, where k-eT denotes the thermal 
energy. This edge tension F is sufficiently large to pre¬ 
vent membrane rupture in this study. Molecular dynam¬ 
ics with a Langevin thermostat is employed42ii^ In addi¬ 
tion to canonical ensemble simulations, replica exchange 
molecular dynamics^i^ for the rod curvature Crod^ is 
used to obtain the thermal equilibrium states. The sim¬ 
ulation data of the replica exchange method (canonical 
simulations) are indicated by lines (the symbols with 
lines) in the figures. The error bars are estimated from 
three or four independent runs. The results are displayed 
with the rod length /rod = lOcr for the length unit and 
k^T for the energy unit. 

The rod assembly on membrane tubes with N = 2400 
and vesicles with N = 2400 and 9600 were investigated 
at the rod density ^rod = N^odNsg/N = 0.1 to 0.5, where 
A^rod is the number of the rods. The tube length is fixed 
as Lz = 4.8rrod in the longitudinal (z) direction. The 
periodic boundary conditions are employed. The radius 
of the tube is i?cyi = 0.989rrod; the radii of the vesicles 
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FIG. 1: Snapshots of a membrane tube with protein rods. 
(a( Grod ~ 0 and 0rod ~ 0.1. (b) GrodCrod ~ 2 and 0rod ~ 

0.1. (c) GrodCrod — 3.2 and (j^rod — 0.1. (d) Grod^rod — 3.75 
and (j^Tod — 0.1. (e) GrodCrod ~ 3.75 and (j^rod ~ 0.25. (f) 
Grod^rod = 3.75 and (p-cod = 0.4. A protein rod is displayed 
as a chain of spheres whose halves are colored in red and in 
yellow. The orientation vector Ui lies along the direction from 
the yellow to red hemispheres. Transparent gray particles 
represent membrane particles. 


are i?ves = 1.54rrod and 3.07ri.od at = 2400 and 9600, 
respectively, in the absence of the rods. 


III. MEMBRANE TUBES 

A membrane tube extended from a vesicle by op¬ 
tical tweezers and micropipette is a typical experi¬ 
mental setup to investigate protein sorting onto lipid 
membranes In our previous paper—, we reported 

the coupling between self-assembly of the protein rods 
and the shape deformation of the membrane tube at low 
rod densities from ((ij-od = 0.1 to 0.25. First, we review 
the rod assembly at low rod density and then describe 
the behavior at higher density. 

At small rod curvatures for 0 < Grod^’rod ^ 2, the 
rods are randomly distributed on the membrane [see Figs. 
[ira),(b)]. As the rod curvature Grod increases, the rod 
orientation changes from the tube (z) axis direction into 
the azimuthal (9) direction, so that the orientational or¬ 
der parameter Sz = /N^od)'^ii'^Si^z'^ — 1) decreases, 
where Si^z is the z component of the orientation vector 
of the f-th rod [see Fig. [2][a)]. With increasing ^rod, 
this orientation change is slightly enhanced by rod-rod 
interactions. 

In a cylindrical tube of a homogeneous membrane, an 
axial force 



is generated by the membrane bending energy, since an 
increase in the axial length results in a decrease in the 
cylindrical radius— Figure [2][b) shows the force fz — fo 
induced by the rods with normalization by (j)rod-, where 
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FIG. 2: Rod curvature Crod dependence of (a) the mean 
orientation degree Sz and (b) the mean axial force fz of the 
membrane tube for = 0.1, 0.167, 0.25, and 0.4. The 
dashed line is obtained from Eq. with Ki = dOfesT. Error 
bars are displayed at several data points. 


/o = 102 kBT/rrod is the force in the absence of the rods. 
All curves in Fig. [2ljb) overlap for 0 < C'rod^’rod ^ 2. 
Thus, the effects of the rods on the axial force are pro¬ 
portional to the rod density 4 >rod unless the rods begin 
to self-assemble. During the rod orientation change for 
0 < C'rod?'rod ^ 1.2, the force fz is almost constant, 
so that it does not contribute to the effective sponta¬ 
neous curvature of the membrane. After the orientation 
is almost saturated (1.2 < C'rod^'rod ^ 2), fz linearly de¬ 
creases with increasing Crod- We will discuss the slope of 
this decrease later in Sec. IlYl 

As Crod increases further (2 < Crod?'rod ^ 3), the rods 
assemble in the azimuthal ( 0 ) direction and the tube de¬ 
forms into an elliptical shape [see Fig. HKc)]. As Crod 
increases even further, the rods also assemble in the longi¬ 
tudinal (z) direction [see Fig. [I][d)]. This two-step assem¬ 
bly is captured by the amplitudes of the Fourier modes 
(see Fig. ED. The lowest Fourier modes of the membrane 
shape and rod density along the longitudinal (z) direc¬ 
tion are given by r^z = (1/lV) ’’’i exp(—27rZii/Tz) and 

nqz = (l/Wod) exp(-27rzii/Tz), respectively, where 
''’f = + Vi- la the azimuthal (0) direction, the m- 

th Fourier modes are Tqgrn = (1/A^) G exp(—mSp) 
and nqSm = (1/Wod) Si exp(-m6»ii), where 0* = 

tada~^{xi/yi). The modes of m = 2 and to = 3 indicate 
elliptical and triangular deformations, respectively. The 
amplitudes of the membrane shape rqe 2 and rod density 
nq 02 along the 6 direction increase together, and subse¬ 
quently, the amplitudes of rqz and Uqz along the z direc¬ 
tion increase at 0rod =0.1. 

As (ftod increases, azimuthal assembly occurs at 



FIG. 3: Rod curvature Crod dependence of the mean ampli¬ 
tudes of (a), (c) shape deformation and (b), (d) rod densities 
of membrane tubes for (^rod = 0.1, 0.167, 0.25, and 0.4. The 
amplitudes of the lowest Fourier mode along the azimuthal 
(6) and longitudinal (z) directions are calculated for the mem¬ 
brane shapes (rqe 2 and r-qz) and densities {nqe 2 and riqz) of 
the center of mass of the protein rods. The Fourier ampli¬ 
tudes are normalized by those at Crod = 0 (denoted by *). 
Error bars are displayed at several data points. 

smaller values of Crod [see Figs. |3][a), (b)]. In contrast, 
longitudinal assembly is suppressed at 4>rod ^ 0.3. This 
is caused by too high density of the rods at the disk edges 
of the tubes [see Figs.[T][e), (f)]. At ^rod = 0.4, the edges 
are completely filled by the rods. The small values of 
< \n*z;\'^ > indicate that little space is left for the rods to 
move translationally along the edges. 

As 0rod increases further, the phase behavior drasti¬ 
cally changes. At tpmd = 0.5, a triangular tube is formed 
with increasing Crod (see Figs. |3] and®). With further in- 
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FIG. 4: Snapshots of membrane tubes at = 0.5. (a) 
Triangular tube at Crod^rod = 3. (b) Triangular tube at 

Grod^rod = 3.4. One of the edges of the triangle is partially 
broken and a discoidal bud is formed, (c) Buckled discoidal 
tube at Grod^rod = 3.25. 


creasing Grod, the triangular shapes become unstable so 
that a buckled discoidal tube is formed [see Fig. |4l[c)]. 
At the largest values of Grod for the triangular tubes 
(Grod^'rod — 3.4), oiie of the edges is not connected by 
the periodic boundary and the rods assemble into a dis¬ 
coidal bud [see Fig. |4l[b)]. 

The resultant tube shapes have hysteresis. With de¬ 
creasing Grod, the buckled discoidal tube can remain at 
Grod^’rod ^ 2.85. Thus, the first-order transition between 
the triangular and discoidal tubes occurs and the tran¬ 
sition point is approximately GrodJ'rod — 3. Since the 
replica exchange method cannot be applied to discrete 
transitions with a large free-energy gap, the exact tran¬ 
sition point cannot be determined in the present simula¬ 
tion. The triangular and discoidal tubes exhibit large am¬ 
plitudes of and respectively [see Figs.[5][a), 

(b)]. Large amplitudes of indicate buckling [see 

Fig.EIc)]. 


IV. SIMPLE GEOMETRIC MODEL OF 
POLYGONAL TUBES 

To understand the dependence of the axial force /z 
and the formation of polygonal tubes, we analyze the 
bending energy Fcv of the membrane tubes in simplified 
geometries. The bending rigidity of the rods parallel and 
perpendicular to the rod axis are denoted as ki and K 2 , 
respectively. The bending rigidity of the membrane in 
the absence of the rods is kq = ISfceT. When the rods 
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FIG. 5; Fourier amplitudes of shape deformation of mem¬ 
brane tubes. The solid lines and symbols represent the results 
of the replica exchange simulation and canonical simulations, 
respectively, at 0rod = 0.5. The discoidal and triangular tubes 
are indicated by ’dis’ and ’tri’, respectively. The dashed lines 
represent the results of the replica exchange simulation at 
0rod = 0.4. Error bars are displayed at several data points for 
the replica exchange simulations. 

are uniformly distributed in a cylindrical membrane tube 
and completely aligned in the azimuthal direction, the 
bending energy Fcv is expressed as 

Fc. = 

( 2 ) 

where A = 2'KRcy\Lz is the membrane area. The axial 
force /z = (9Fcv/i9LzU is given by 

fz = + ‘^^4>Tod “ S^lC'rod) • (3) 

The first term is the force (/o) in the absence of the rods. 
The second term is the force generated by the rods, which 
is proportional to 4>rod and linear with respect to Grod- 
These dependences agree with the simulation results (see 
Fig. [2]). Note that one should not take preaverages of 
the bending rigidity and spontaneous curvature in Fcv 
(k = Ko(l - 4>Tod) + Ki(j)rod and Gq = Grod^rod)- If they 
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are preaveraged, the effects of the rods are overestimated 
and fz is not a linear function of ^rod ■ The bending rigid¬ 
ity Ki is estimated as ki = AOk^T = 2.67kq from the fit 
of the slope of /^ (see the dashed line in Fig. [5)). This 
value is reasonable for the parameter set of our simula¬ 
tion. If hexagonal packing is assumed around an isolated 
protein rod, the protein particle has four neighboring 
membrane particles and two protein particles. Although 
the coefficient of the bending potential along the rod is 
four times larger than that between membrane particles, 
the effective bending rigidity is slightly reduced by the 
membrane-rod interactions parallel to the rod. The ab¬ 
solute value of /^ obtained by Eq. @ is smaller than the 
simulation results. This is likely caused by the assump¬ 
tion of the complete azimuthal alignment of the rods and 
the negligence of the weakly attractive interactions be¬ 
tween the rods. 

Next, we consider simplified geometries of polygonal 
tubes as shown in Fig. [61(a). The n-gonal tube consists of 
n flat membranes [black lines in Fig. IHKa)] and n round 
edges with a length Trod and a curvature 27r/nrrod [red 
curves in Fig. EKa)]. All of the protein rods are assumed 
to be aligned in the azimuthal direction and uniformly 
distributed in the edge regions. The n-gonal tube can 
exist as a straight tube at A(j).[od < Aed < A, where 
the edge area Aed = nr^odAz and the membrane area 
A = 2TTRcy\Lz. The edges are buckled at At/frod > Agd- 
The membrane does not have sufficient area to form an n- 
gonal tube at Agd > A (i?cyi/?'rod < n/27r). The bending 
energy of the straight polygonal tube is given by 

Ecvl = -rr (-C'rod ) A())rod (4) 

2 V nr rod 7 

2 Vnrrod^ 

The buckled polygonal tube also has non-zero cur¬ 
vature along the longitudinal direction. Here, we as¬ 
sume that the curvatures of the edge region are constant 
and the curvature energy in the other principal direction 
(close to the azimuth) is approximated as Eq. (|1|) with 
Agd = A())rod = ncrodfred, where Led is the edge length 
and Led > L^. The curvature energy along the longitu¬ 
dinal direction is given by 


Fcv2 — 


32/c^nrrod ^ E{k)K{k)] (5) 

32K2nrrod f 2 

'ly 




T I S ■ ■ Rd - (6) 

Lgd I 8 64 J 

dTT^Kanrrod r _ 1^^ _ _^^2'1 

1 Led 4 Led 1 


Led 


where K{k) and E{k) are the complete elliptic integral of 
the first kind and the second kind, respectivelyi^^i^ All 
of the edges are buckled with the same amplitude. The 
modulus k is determined from the length ratio. 


^ 2E{k) ^ 

Led K{k) 

= 1 - - fc4/8-ko(fc®). 


( 8 ) 



FIG. 6: Membrane tube structures obtained by a simple 
geometric model, (a) Schematic representation of slices of 
polygonal tubes with the number of edges n = 2, 3, and 
4. (b) Curvature free energy Fcv of polygonal tubes at 

Acyi/rrod = 1, (/>rod = 0.5, aud Lz/r^od = 4.8. (c) Phase 

diagram at Rcyi/r^od = 1- The solid lines represent the phase 
boundary at Lzlr^od = 4.8. The dashed line represents the 
phase boundary between the discoidal tube (n = 2) and tri¬ 
angular tube (n = 3) at Tz/crod = 2.4. 


For small buckling amplitudes, F'cv 2 can be expanded for 
1 - Tz/Lgd <C 1 as in Eq. d?]). 

Figure EKb) shows an example of Fey obtained by 
Eqs. ([3]) and ([7]) with ki = 2.67ko and K 2 = kq. At 
larger values of Crod, the polygonal tube with a smaller 
n has the lowest energy. The phase diagrams for (p^od 
versus C^d at the simulation condition are shown in 
Fig.EKc). The phase boundary is slightly underestimated 
(C'rodCrod = 2.6 in the simple model and Crod?'rod — 3.2 
in the simulation). This deviation is caused by neglecting 
the translational and orientational entropies of the rods. 

The polygonal edges are buckled at prod > 
nrrod/‘2'!TRcyi. At the buckling transition point, the phase 
boundary exhibits a kink. The buckling of a discoidal 
tube with n = 2 starts at prod = 0.32. The phase bound¬ 
ary between the n = 2 and n = 3 phases increases as 
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FIG. 7: Phase diagrams of polygonal tubes obtained by 
a simple geometric model for (a) Crod^rod = 3 and (b) 
Grod^rod — 1-9 at Lzj'^YoA — 4.8. 


the buckling amplitude increases with increasing 4>Yod [see 
Fig.EIc)]. As the triangular tube with n = 3 is buckled 
at d'rod > 0.48, the boundary decreases again. With a 
decrease in the tube length Lz, the buckling amplitude 
increases so that the phase boundary exhibits a greater 
increase, whereas the boundary between straight polygo¬ 
nal tubes has no Lz dependence [compare the dashed and 
solid lines in Fig. [DJc)]. Reentrant phase transitions oc¬ 
cur when the boundary of buckled and straight polygonal 
tubes is crossed with increasing ^rod- 

In experiments, protein density and tube radius are 
controllable parameters. Figure [7] shows the phase dia¬ 
gram for (ji-Yod versus i?rod at fixed values of Grod- As the 
tube radius increases, tubes with large values of n can 
exist. At large curvatures of Grod, the transitions occur 
between the straight tubes at small rod densities. In con¬ 
trast, transitions between buckled polygonal tubes can be 
more complicated at 1 < Grod ^2.5. At Grod?’rod = 1.9, 
reentrant phase transitions occur between n = 3 and 4 
[see Fig. |7l[b)]. With increasing a buckled triangu¬ 
lar tube transforms into a straight square tube. Subse¬ 
quently, the square tube is buckled and transforms back 
to a buckled triangular tube. 



FIG. 8; Vesicle shapes at c^rod = 0.167 for N = 2400 and 
9600. The snapshots are shown at (a) Grod^rod = 1.5, (b) 
2.5, and (c) 3.5 for N = 9600. Top and bottom snapshots 
are displayed in bird’s-eye and front views, respectively, (d) 
Three mean eigenvalues of the gyration tensor of the vesicle. 
The dashed and solid lines represent data for N = 2400 and 
9600, respectively. 


V. VESICLES 

Self-assembly of the protein rods also occurs in a vesi¬ 
cle. At small values of Grod, the rods are randomly 
distributed on the vesicle [see Fig. |Sl[a)]. With an in¬ 
crease in Grod, the rods assemble into an equator of the 
vesicle (parallel to the rod axis) and the vesicle deforms 
into an oblate shape [see Fig. [El[b)]. With a further in¬ 
crease in Grod, the rods also assemble along the equator 
(perpendicular to the rod axis). Thus, the rods assem¬ 
ble in both directions, and a discoidal bump is formed 
[see Fig. HKc)]. In our previous paper—, we reported the 
rod assembly at (ji^od = 0.167 and N = 2400. Figure [8| 
shows that this assembly also occurs with large vesicles 
at TV = 9600. These shape deformations are quantified 
by the changes of three eigenvalues of the gyration ten¬ 
sor [see Fig. HKd)]. The gyration tensor is expressed as 
a-ap = -aG){l3j -/3g), where a,l3€x,y,z 

and ac is the center of mass. When the vesicle is ap¬ 
proximated as an ellipsoid, these eigenvalues are squares 
of the three principal lengths of the ellipsoid. For the 
deformation into the oblate shape, the smallest length z 
decreases and the other lengths x and y increase. For the 
bump formation, the largest length x increases. 

At high rod density 4>Yod, it is found that polygonal 
vesicles are formed in a manner similar to that of the 
polygonal tubes (see Fig. o. As 0rod increases, the vesi¬ 
cle equator is completed filled by the rods so that phase 
separation along the equator does not occur. With a fur- 
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FIG. 9: Vesicle shapes for various (f>rod values a.t N = 9600. 
(a) Snapshots of vesicles of an elliptical disk and a triangular 
hosohedron at Crod^'rod = 2.5 and 0rod = 0.3. (b) Snap¬ 

shots of vesicles of a tetrahedron and a triangular prism at 
Grod^rod = 2.5 aud (piod ~ 0.4. (c), (d) Three mean eigenval¬ 
ues of the gyration tensor of the vesicle. The labels ’dis’ and 
’tri’ denote the vesicles with disk and triangular hosohedral 
shapes, respectively, (c) Dependence on Grod at </)rod = 0.3. 
(d) Dependence on (jj-^od at GrodDod = 2.5. 


ther increase of (^rod, the vesicle deforms into an elliptic 
disk instead of buckling [see the top snapshot in Fig.jSlJa) 
and the increase of in Fig. [5][d) from (/j^od = 0.25 to 
0.35]. With an even further increase, the vesicle trans¬ 
forms into a triangular hosohedron [see the bottom snap¬ 
shot in Fig. [5][a)]. The elliptic disk and triangular hoso¬ 
hedron coexist for finite ranges of i^rod and Crod [see 


Figs.lUd) and (c)]. These are the same types of behavior 
between discoidal and triangular tubes. For (pj-od = 0.4, 
the tetrahedron and triangular prism are obtained as 
metastable states. The energy analysis in Sec. |TV] sug¬ 
gests that these shapes and polyhedra with more edges 
can be stabilized at larger vesicles. 

For polyhedra with genus 5 = 0, the number of 
faces rif, vertices riv, and edges rie have the relation 
n.f + Uv — u-e = 2 (Euler’s formula). In the polyhedra 
induced by the rods, three edges are connected at any 
vertex. This is different from elastic capsules^^— that 
prefer an icosahedron, where five edges are connected at 
the vertex. Since each edge contacts two vertices, the 
relations Ue = and rie = 3(nf — 2) are obtained. 

Fluid vesicles with linear defects^ also have this type of 
polyhedra. 


VI. SUMMARY 

We have revealed that the formation of polygonal tubes 
and polyhedral vesicles can be induced by banana-shaped 
protein rods. At low rod density, rod self-assembly occurs 
in two steps. As the rod curvature increases, the rods as¬ 
semble in the direction parallel to the rod axis and form 
discoidal tubes and vesicles. Subsequently, the rods as¬ 
semble along the disk edges. As rod density increases, the 
phase behavior remarkably changes. With an increase in 
the rod curvature, the membrane tube transforms into 
a triangular tube and subsequently form a buckled dis¬ 
coidal tube. The vesicle transforms into an elliptic disk 
and triangular hosohedral shapes. The transitions be¬ 
tween triangular and discoidal shapes are the first-order 
phase transitions. These transitions can be understood 
by simple geometric analysis of the bending energy. 

In our simulations, the number of protein on the mem¬ 
brane is fixed. In typical experimental condition, the 
proteins are inserted from the bulk so that the protein 
density on the membrane is determined by the balance of 
the chemical potentials. Thus, in experiments, the pro¬ 
tein density may discretely change accompanied by the 
shape transition between the polygonal tubes (or pohy- 
hedral vesicles). The polygonal tubes with more faces 
can be adsorbed by more proteins. The drastic change of 
the the protein density may indicate a shape transition. 

Polygonal tubes and polyhedral membrane shapes are 
often observed experimentally. Regular arrays of triangu¬ 
lar prismatic tubes are observed in the inner membranes 
of the mitochondria in astrocytesi^^— The formation 
mechanism of these triangular tubes is not known. Our 
simulation study suggests that these prismatic tubes can 
be generated by the assembly of the BAR proteins. 

In the presented simulation conditions, tubulation 
is not obtained. Tubulation is likely prevented by the 
high bending energy at the tip of the membrane tubule, 
resulting in linear rod aggregation at the edges of the 
polyhedra. The addition of the spontaneous curvature 
of the membrane induces tubule formationi^ In living 

















cells, the tubules and discoidal membranes are often 
connected. Endoplasmic reticulum (ER) consists of flat 
discoidal membranes (rough ER) with branched tubular 
network (smooth ER). It is an interesting problem to 
clarify the control mechanism of tubules and discoidal 
membranes for further studies. 
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